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Abstract. In this paper, we recalibrate the Cepheid dis- 
tance to some nearby galaxies observed by the HST Key 
Project and the Sandage-Tammann-Saha group. We use 
much of the Key Project methodology in our analysis but 
apply new techniques, based on Fourier methods to esti- 
mate the mean of a sparsely sampled Cepheid light curve, 
to published extra-galactic Cepheid data. We also ap- 
ply different calibrating PL relations to estimate Cepheid 
distances, and investigate the sensitivity of the distance 
moduli to the adopted calibrating PL relation. We re- 
determine the OGLE LMC PL relations using a more 
conservative approach and also study the effect of using 
Galactic PL relations on the distance scale. 

For the Key Project galaxies after accounting for 
charge transfer effects, we find good agreement with an 
average discrepancy of -0.002 and 0.075mag. when us- 
ing the LMC and Galaxy, respectively, as a calibrat- 
ing PL relation. For NGC 4258 which has a geometric 
distance of 29.28ma<?., we find a distance modulus of 
29.44 ± 0.06(random)mag., after correcting for metallic- 
ity. In addition we have calculated the Cepheid distance 
to 8 galaxies observed by the Sandage-Tammann-Saha 
group and find shorter distance moduli by — 0.178mag. 
(mainly due to the use of different LMC PL relations) and 
— 0.108map. on average again when using the LMC and 
Galaxy, respectively, as a calibrating PL relation. However 
care must be taken to extrapolate these changed distances 
to changes in the resulting values of the Hubble constant 
because STS also use distances to NGC 3368 and 4414 
and because STS calibration of SN la is often decoupled 
from the distance to the host galaxy through their use 
of differential extinction arguments. We also calculate the 
distance to all these galaxies using PL relations at maxi- 
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mum light and find very good agreement with mean light 
PL distances. 

However, after correcting for metallicity effects, the 
difference between the distance moduli obtained using the 
two sets of calibrating PL relations becomes negligible. 
This suggests that Cepheids in the LMC and Galaxy do 
follow different PL relations and constrains the sign for 
the coefficient of the metallicity correction, 7, to be neg- 
ative, at least at the median period log(P) w 1.4, of the 
target galaxies. 

Key words: Cepheids - distance scale 



1. Introduction 

The extra-galactic distance scale is one of the key prob- 
lems in modern astronomy. One of the basic parts of the 
solution is the correlation between period and mean lu- 
minosity (PL) obeyed by classical Cepheids. The Hubble 
Space Telescope (HST) H Key Project (Freedman et al. 
2001, hereafter KP) has used Cepheids (discovered either 
by ground-based telescopes or by HST) in 31 spiral galax- 
ies, with 18 of them observed originally by KP, and the PL 
relation to estimate the distances to these galaxies. These 
distances were then used in turn to calibrate a host of sec- 
ondary distance indicators and hence estimate Hubble's 
constant. In parallel with this a team led by Sandage et 
al. (Saha et al. 2001b, hereafter STS) has used the HST 
to discover Cepheids in spiral galaxies which host Type la 
supernovae. Cepheid distances to these galaxies were used 
to calibrate the Hubble diagram for Type la supernovae 
and hence estimate Hubble's constant. 

Though both groups used mostly HST observations to 
discover the Cepheids and similar methodologies, the KP 
team favor a short distance scale and a larger value of 
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H (72±8kms- 1 Mpc- 1 , Freedman et al. lMiT)) whilst the 
STS group favor a long distance scale and smaller value of 
H (58.7±6Mms- 1 Mpc-\ Saha et al. l2001bll . The dis- 
crepancy in the value of the Hubble constant from these 
two groups is still unresolved (e.g., see Hendry & Kanbur 
IT5551 Beaulieu et al. HMTl Kochanek HMTI Sasselov et al. 
HM7I Tanvir et al. lT555l Caputo et al. fMM Gibson et al. 
I2UUU1 Gibson & Stetson l2Uo"T)l . It is due only in part to the 
Cepheid distance scale. We note that the recent WMAP 
results (Spergel et al. find a Hubble constant very 

similar to the value obtained by KP. Nevertheless it is still 
very instructive to discover the reason behind the discrep- 
ancies in the Cepheid distance scale from the STS and 
KP groups. A proper understanding of the Cepheid PL 
relation in this regard is very important for an accurate 
local distance scale. In this paper we concentrate on re- 
calibrating the Cepheid distance to the target galaxies in 
both KP and STS groups with existing data because of 
the following factors: 

1. HST Cepheid data is sparsely sampled with typically 
12 and 4 points per Cepheid in the V and I band re- 
spectively (see Section 2 for details) . In order to apply 
the Cepheid PL relations, it is necessary to estimate 
the mean magnitudes from these data in both V and I 
bands. In this paper, we use the techniques of Fourier 
expansion and interrelations to estimate the V and I 
mean magnitudes from sparsely populated light curves. 
These techniques arc described in detail in Ngeow et 
al. plol)) . 

2. We re-analyze the OGLE (Optical Gravitational Lens- 
ing Experiment) LMC Cepheids (Udalski et al.|l999a) 
on the basis of the quality of their V band light 
curve, and develop more conservative LMC PL rela- 
tions based on this analysis and investigate the sensi- 
tivity of extra-galactic Cepheid distances to the LMC 
PL relation. This approach is different to the "sigma- 
clipping" methods currently used by Udalski et al. 

mm- 

3. The average value of metallicities, defined as 12 + 
log(0/if), in target galaxies is about 8.84 ± 0.31dex 
(Freedman et al. 12001(1 . This value is closer to the stan- 
dard Solar value of 8.87±0.07dex (Grevesse et al. H55fy 
than the LMC value of 8.50 ± 0.08dex (see reference 
in Ferrarese et al. 12000(1 . Therefore, another approach 
is to use Galactic Cepheid PL relations as fundamen- 
tal calibrating relations (Feast 2001] 120031 Tammann 
I2UU51 Thim et al. 1200^1 Fouque et al. 11003) ■ We com- 
pare the distances obtained when using both the LMC 
and Galactic PL relations. 

4. Kanbur & Hendry (1996) building on the hydrodynam- 
ical models of Simon et al. (1993) suggested that PL 
relations at maximum light may be more accurate than 
PL(mean) relations. To test this idea, we compute dis- 
tances using PL (max) relations and compare with their 
mean light counterparts. 



In Section 2 we summarize the photometric data used 
in this study. In Section 3 we present and discuss the 
methodology used in determining Cepheid distances, in- 
cluding obtaining the V and I band means and applying 
different PL relations. The results will be presented in Sec- 
tion 4. The conclusion and the discussion will be presented 
in Section 5. 

2. The Photometric Data 

The photometry data for the Cepheids in each target 
galaxy were directly obtained from the corresponding pub- 
lished papers. We emphasize that we did not repeat the 
photometric reduction from raw data. The target galaxies 
that are selected in this study include: 16 KP galaxies, 8 
STS galaxies and NGC 4258, which has an accurate geo- 
metrical distance measurement from water maser studies 
(Herrnstein et al. 119991 hereafter WM galaxy) . The pho- 
tometric data of 16 KP galaxies were directly downloaded 
from the KP web-page 1 , excluding NGC 1425 (the data 
is not available at the time of analysis) and NGC 5457 
(M101, as the observations to this galaxy include its outer 
field (Kelson et al. 11996(1 and inner field (Stetson et al. 
1998b,), which complicated the analysis). The photomet- 
ric data for the STS and WM galaxies were taken from the 
STS papers (see Saha et al. I2001bl for references for each 
target galaxy) 2 and Newman et al. 1(2001(1 . respectively. 
The list of the target galaxies can be obtained from Table 

There are two major photometry reduction pack- 
ages used in reducing the data for these galaxies, the 
ALLFRAME (Stetson IT5M1 UMB) and a variant of the 
DoPHOT (Schechter et al. IT5551 Saha et al. HMfy pack- 
age. The KP team utilized both packages as a double- 
blind reduction process (Kennicutt et al. 119951 Freedman 
et al. 12001(1 to check the consistency between the two 
packages. However, the final, adopted distance moduli are 
based on ALLFRAME results (Freedman et al. l2"Uo"Tj) . 
Inter-comparison of the results from these two packages 
(galaxy-by-galaxy basis) show good agreement between 
the two (see the KP papers for more details and Hill et 
al. 11555)1 . We also use the ALLFRAME photometric data 
downloaded from their web-page (if available). However, 
the photometric data for NGC 2541 and NGC 4321 are 
only available in DoPHOT from the same web-page. Sim- 
ilarly, only the DoPHOT photometric data is available 
for WM galaxy, although the photometry reductions were 
carried out by both packages (Newman et al. 2001). 

In contrast, the photometric data for the STS galax- 
ies were reduced mainly with DoPHOT, although some of 

1 The URL is http://www.ipac.caltech.edu/H0kp/H0KeyProj 
.html. The reference to each KP galaxies also listed in this 
URL. 

2 We exclude NGC 3368 (Tanvir et al. 115551 because it's not 
part of the STS program, and we are not calibrating the SN la 
distance in this paper. 
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the galaxies also used an additional reduction package to 
check the DoPHOT results (e.g., reduction of NGC 4527 
included the ROMAFOT package, see Saha et al. l2001a|) . 
Nevertheless, the available data for all STS galaxies are 
from the DoPHOT package. Although the KP team has 
reanalyzed the STS galaxies with the ALLFRAME pack- 
age (Gibson et al. ESTOH Gibson & Stetson IMiTl Stetson 
& Gibson 2001] Freedman et al. I2001[l . these data are not 
yet available from the KP web-page (as claimed in these 
papers). However, most, but not all, of the Cepheids are 
common in both STS and KP results (Gibson et al. 2000). 
Hence, we only can use the DoPHOT results in our anal- 
ysis of STS galaxies. 

The next steps in the photometric reduction process 
are a Cepheid search and period determination. Each re- 
duction package has its own algorithms to perform the 
period search, and the results generally agree well. We 
use published Cepheid data and periods in our analysis. 

The HST observations of the target galaxies generally 
consist of 24 V band and 8 I band cosmic-ray (CR) split 
images. The analysis of DoPHOT and ALLFRAME treat 
these images differently to deal with the cosmic-ray and 
point-spread functions. However, most of the published 
photometry combined the back-to-back CR split images 
together and tabulated the averaged 12 V band and 4 I 
band data points. Therefore, we only can use these (re- 
duced) photometric data points to reconstruct the light 
curves and obtain the mean magnitudes by the Fourier 
techniques described in the next section. 

3. General Methodology and Analysis 

The physical basis of Cepheid PL relations and their us- 
age in determining the distance have been covered exten- 
sively in the literature (e.g., see Sandage & Tammann 
HMHI Feast k Walker H%71 Madore & Freedman ITMT1 
Freedman et al. 2001) and would not be repeated here, 
only a brief description will be presented. The application 
of the Cepheid PL relation to estimate distances involves 
the discovery and appropriate observation of Cepheids suf- 
ficient to estimate their periods (in days) and mean mag- 
nitudes. In order to correct for the extinction/reddening, 
observations of extra-galactic Cepheids by HST are nor- 
mally taken in bandpasses V and I. The PL relation in 
bandpass A can be expressed as: M\ = a\ log(P) + b\ 7 
where a and b are the slope and zero-point, respectively. 
For the measurement of the apparent mean magnitudes of 
Cepheids in a target galaxy, < m\ >, the distance modu- 
lus in bandpass A is: 

Ma =< m x > -a A log(P) - b\. 

Since /io = /iy — Ay = fii — Aj, the reddening-free dis- 
tance modulus, or the Wessenheit function (Madore 1982 
Moffett & BarnesHHEHl Madore & Freedman HTMl Tanvir 
1999]), can be derived as: 



fi = Mv - - A*/), (!) 

where R = Av/E(V — I) is the ratio of total-to- 
selective absorption. Following KP, we adopt R = 2.45 
from Cardelli et al. fl989j) . Note that the validity of equa- 
tion (1) is based on the assumption that there are no corre- 
lated measurement errors in the V and I bandpasses. Then 
any differences in the V and I band distance moduli are 
due to differential reddening (see Saha et al. 119961 12001bl 
for details). 

After applying a period cut to the short period 
Cepheids in target galaxies (to avoid the incomplete- 
ness bias at the faint-end of observations, see Sandage 
[T^HHI Lanoix et al. 11999a! Freedman et al. |2"UUT|) . an un- 
weighted mean of apparent distance moduli to the remain- 
ing Cepheids is adopted as the distance modulus to the 
target galaxy. The random (or statistical) error associated 
with equation (1) can be calculated via the standard for- 
mulae (i.e. a/y/N). The reasons for taking the unweighted 
mean are: (a) the photometric errors of the mean magni- 
tudes are smaller than the expected width of the insta- 
bility strip, and are, together with other systematic er- 
rors, the dominant part of the weight in a weighted mean 
(Leonard et al. [5D03 ,) ; (b) this is equivalent to the fitting 
scheme used by KP; and (c) it can be easily incorporated 
with the weighting scheme adopted by the STS group 
fTanvir H997[) . Finally, due to the possible metallicity de- 
pendence of the Cepheid PL relation (see the discussions 
and reference in, e.g., Kennicutt et al. 119981 Freedman 
et al. I2001fl . metallicity corrections are added to equation 
(1): 

Hz = Ho + S z , (2) 

where S z = j([0/H] re f — [0/H] ga i) with the usual def- 
inition of [O/H] = 12 + \og{0/H). The reference metal- 
licity, [0/H] re f, is 8.50dex and 8.87dex for using LMC 
and Galactic PL relations, respectively. Some published 
values of 7 range from -0.88 ± 0.16 and -0.56 ± 0.20 
(Goul d Pip , — Q.44±g-|g (Beaulieu et al. ITM71 Sasselov 
et al. IT39T|) . -0.4 ± 0.2 CKochanek HMTjl . -0.32 ± 0.21 
(Freedman & Madore lTMll . -0.24±0.16 (Kennicutt et al. 
HM^ to 0.27 (Caputo et alEES- Note that the value of 
7 is method-dependent, i.e. it depends on the bandpasses 
and calibrating PL relation adopted. For example, strictly 
speaking the metallicity dependence quoted in Kennicutt 
et al. l|1998fl is only applicable when the Madore & Freed- 
man Ijl991|l PL relation is used, and then only for the V 
and I bandpasses. Although there is some debate on the 
value for 7 (see Freedman et al. 120011 Tammann et al.]2001 
for details), we adopt 7 = jv.i = —0.2±0.2mag dex~ x as 
in KP, who note that this value is in the middle of a num- 
ber of different determinations of 7. Changes in metallicity 
affect the mean brightness of a Cepheid and the papers 
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quoted above were aimed at quantifying this brightness 
shift due to mctallicity. However, we now know that the 
Cepheid PL relation has a different slope in the Galaxy 
and LMC, and indeed may have two different slopes in 
the LMC itself. Thus it may be that metallicity affects 
not only the mean brightness of a Cepheid but also the 
way that mean brightness changes with period (the slope 
of the PL relation). We discuss this and its implications 
for metallicity corrections like equation (2) in Section 5. 

Since we are using the published photometry data to 
calculate the distance modulus, we share the same system- 
atic errors as in the published results. These include the 
uncertainty in the zero point of the PL relation, calibra- 
tion of photometric zero points, reddening and metallicity 
effects, incompleteness bias and crowding. The order of a 
typical systematic error is around ~ 0.1 — 0.2mag. 

3.1. The Mean Magnitudes 

As mentioned in Section 2, most of the HST observa- 
tions of extra-galactic Cepheids contain only 12 V band 
epochs and 4 (or 5) I band epochs. These observations 
use a power-law time series sampling strategy to minimize 
the aliasing problem and maximize the phase coverage of 
the observed data points within the observing windows 
(Freedman et al. 119941 Kennicutt et al. I1995|l . In order 
to reduce the bias due to sampling procedures, STS and 
some early KP observations used phase-weighted inten- 
sity means to find mean magnitudes in the V band. The 
I band mean magnitudes were found via empirical rela- 
tions developed by Freedman 1)1988(1 and Labhardt et al. 
(1997) for KP and STS observations, respectively (see, 
e.g., Silbermann et al. 1996). Other KP observations used 
a template- fitting procedure (Stetson I19BSI to obtain the 
mean magnitudes in both V and I band simultaneously. 

In addition to the methods mentioned, Ngeow et al. 
(2003) recently developed an alternative method to ob- 
tain the mean magnitudes, which is based on Fourier tech- 
niques. Ngeow et al. I|2003|) give specific examples of sit- 
uations when such a method could be preferable to tem- 
plate based techniques. Another motivation for developing 
these Fourier techniques is the possibility of reconstruct- 
ing observed HST light curves to compare with models. 
In Section 4.1, we quantitatively compare the means ob- 
tained using our methods with existing techniques in the 
literature. 

Here, we apply the new Fourier method to all target 
galaxies in order to be consistent in doing the analysis, and 
also to test this new method in the distance scale problem. 
The details of using the Fourier techniques are presented 
in Ngeow et al. ((2003 ■ Here we outline the steps: 

1. The V band: The light curves of 12 V band data points 
are reconstructed via a 4 t/l -order Fourier expansion 



with the form: 

i=4 

m(t) = m + y^}Aj cos(2iirt/P + 0»)] (3) 

i=l 

where A and <p are Fourier amplitudes and phases, 
respectively, and too is the mean magnitude. Since 
the periods are directly taken from literature, the re- 
maining nine parameters can be obtained by fitting 
the data with equation (3). We use the technique of 
simulated annealing and restrict the ranges that the 
Fourier parameters in equation (3) can take to recon- 
struct the V band light curves. The ranges of Fourier 
phases are from to 2ir, whilst the ranges of Fourier 
amplitudes are determined from the "calibrating set" 
Cepheids and OGLE LMC Cepheids. The "calibrat- 
ing set" Cepheids consist of ^100 Galactic Cepheids, 
mostly observed originally by Moffett & Barnes ( 1984), 
and some LMC/SMC Cepheids. 

2. The I band: For 4 I band data points, it is clear that 
equation (3) cannot be applied to reconstruct the I 
band light curves. Therefore, we use Fourier interre- 
lations (Ngeow et al. [2003) to reconstruct the I band 
light curves. The Fourier interrelations are linear rela- 
tions connecting the Fourier parameters in the V and 
I bands: 

Ai(I) = ai+PiMV); 4>i(I)=Ji + vMV) (4) 

The coefficients in equation (4) are determined from 
the "calibrating set" Cepheids. Ngeow et al. (2003) 
developed such Fourier interrelations for the Galaxy, 
LMC and SMC separately and found that they are 
only weakly dependent on metallicity. This is impor- 
tant since it means that the interrelations can be ap- 
plied over a wide range of metallicity. Given the so- 
lution of equation (4), the observed I band points are 
then used to establish the I band light curves and hence 
estimate the corresponding mean magnitudes by min- 
imizing the chi-square. 

3. These Fourier techniques have been tested with Monte 
Carlo simulations, and show that the reconstruction 
procedures are unbiased and the errors of the Fourier 
amplitudes and means are around ~ 0.03TOag. 

Although the filters installed on the HST are close 
to the standard Johnson and Kron-Cousins system, i.e. 

~ V and F8UW - I, conversions from the HST 
passbands to the standard photometry systems have been 
used by both KP and STS teams. For the case of the ALL- 
FRAME photometry reduction in the KP galaxies, the 
conversions are made during the reductions with the for- 
mulae suggested by Holtzman et al. ( 1993J), hence the pub- 
lished photometry of the Cepheids are in standard bands. 
A similar situation regarding the conversion holds for the 
two KP galaxies with DoPHOT photometry (NGC 2541 
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Table 1. Mean V and I magnitudes from Fourier fits. 3. 



Galaxy-Cepheid 


log(P) 


V(mag.) 


I (mag.) 


NGC925-C5 


1.686 


23.680 


22.478 


NGC925-C6 


1.635 


24.533 


23.453 


NGC925-C7 


1.624 


24.305 


23.457 



a The entire table is available electronically at the CDS. 



& NGC 4321) and the WM galaxy. For STS galaxies with 
DoPHOT photometric reductions, the conversions are ap- 
plied after the reductions to the calculated mean magni- 
tudes (see equations (5) and (6) in Saha et al. 119941 for 
WFPC; and equations (2) and (3) in Saha et al. HMrjl 
for WFPC2). In the case of PL(max) relations, however, 
no such conversion is available yet. We have to assume the 
conversions at maximum light are similar to the equations 
given by Saha et al. I|19941 [T^fj)l . and apply the conver- 
sion accordingly. These conversions have been applied to 
STS galaxies in Table El 

The derived mean magnitudes in both V and I bands 
for the Cepheids used in this study are listed in Tabled 
where column 1 shows the Cepheid in the target galaxy, 
column 2 is the period adopted from the literature, and 
columns 3 and 4 are the mean V and I band magnitudes, 
respectively, derived from the Fourier techniques described 
in this paper. Tablets available in its complete electronic 
form at the CDS 3 . Here we only show a portion of the table 
to indicate its form and content. 

3.2. The Period-Luminosity Relations 

In order to apply the PL relations to HST data, we only 
look for the published PL relations that are available in 
both V and I bands in both the Galaxy and LMC. The 
following subsections discuss the adopted PL relations in 
this study. Note that for the LMC PL relations, we adopt 
t^LMC = 18.50mag., to be consistent with the KP team 
(Freedman et al. l2U(Jlf) . All of the adopted PL relations are 
listed in Table [3] and [2 for the V and I bands, respectively. 

3.2.1. The LMC PL Relations 

The LMC PL relations used by both KP team (before the 
publication of their final paper) and STS team are based 
on a homogeneous sample of 32 Cepheids, with periods 
ranging from 10 days to ^120 days (Madore & Freedman 
1991, hereafter MF91). Since this PL relation has been 
extensively used in determining Cepheid distances in the 
past, and the STS team still use this PL relation in their 
study (e.g., see Saha et al. 12001 bf) . we adopt the MF91 re- 
lations as one of our calibrating PL relations in this study. 

After the publication of PL relations derived from 
OGLE LMC Cepheids (Udalski et al. Il999al hereafter 

3 http:/ /cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/ 



U99), the KP team recalibrated their Cepheid distances 
with these new LMC PL relations in their final paper, as 
well as the Cepheid distance to NGC 4258 by Newman et 
al. H2001fl . These reddening-corrected PL relations were 
derived using Cepheids with logP > 0.4 to minimize pos- 
sible contamination by first overtone pulsators and were 
" sigma clipped" , resulting in ~ 650 Cepheids used in de- 
riving the PL relations. In order to compare our results 
with published Cepheid distances, we also adopt the U99 
PL relations. The new U99 PL relations have dramatically 
changed Cepheid distances compared to distances derived 
from MF91 PL relations, in the sense that the derived dis- 
tances are smaller with U99 PL relations (Freedman et al. 
1277011 . 

There are some criticisms about the U99 PL relations 
in the literature. The U99 PL relations are dominated by 
short period Cepheids (with < log(P) >^ 0.5 and about 
90% of them have period shorter than 10 days), and a 
lack of longer period Cepheids with log(P) > 1.5 (Feast 
l2T)0H ISnnTn Saha et al. I2001b|l . However, Freedman et al. 
i|2T)0Hl . by using Sebo et al. LMC data, claimed 

that the use of the OGLE LMC PL relations in estimat- 
ing distances to target galaxies whose Cepheids all had 
periods longer than the longest period Cepheid observed 
by Udalski et al. I|1999af) made little difference. Another 
potential problem of the U99 PL relations is the discovery 
of a break in the PL relation at log(P) = 1.0, as shown in 
Tammann et al. ('2001) and in Tammann & Reindl H2002[l . 
The long and short period Cepheids follow different PL 
relations. The reasons for the break in PL relations are 
still unclear, nevertheless we adopt the PL relations for 
long period Cepheids (Tammann & Reindl 2002 , hereafter 
TR02) since most of the extra-galactic Cepheids have pe- 
riods longer than 10 days. 

Despite these potential problems, we re-analyze OGLE 
LMC Cepheid data and re-fit the PL relations without us- 
ing sigma-clipping. The sigma-clipping method is an iter- 
ative procedure whereby the data are fitted with a regres- 
sion and then those points lying 2.5cr away from the regres- 
sion line are removed. A new regression is fitted and the 
procedure continues for few cycles. Therefore, it blindly 
removes all the outliers to reduce the scatter in the fi- 
nal fit. However, we feel this approach may both remove 
some true Cepheids from the data and include some suspi- 
cious Cepheids in the final sample (see discussion below). 
In addition, Nikolaev et al. I|2()()3p discussed some flaws 
associated with the sigma-clipping algorithm, including 
the implicit assumption of the normal distribution of the 
residuals and the sensitivity of the results to the chosen 
threshold of ka. 

We use the same OGLE LMC data as in Udalski et 
al. 1|1999a[) to derive the LMC PL relations, but with- 
out processing the sigma-clipping procedure (Udalski et 
al. 11999a! Willick & Batra l2TjTT01 Groenewegen l27j7J0|) . 
Instead, we examine in detail the light curves for all 
Cepheids. The OGLE photometric data were downloaded 
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Table 2. The adopted V band Period-Luminosity relations. a 



Relation 


Slope (fly) 


ZP (b v ) 


av 


N 


LMC (MF91) 
OGLE LMC (U99) 
OGLE LMC (Here) 
LMC PL >10d (TR02) 
LMC PL (Max) 


-2.76 ±0.11 
-2.760 ±0.031 
-2.746 ± 0.043 

-2.48 ±0.17 
-2.744 ±0.051 


— 1 40 ± 05 
-1.458 ±0.021 
-1.401 ± 0.030 

-1.75 ±0.20 
-1.817 ±0.035 


0.27 
0.159 
0.223 
0.16 
0.261 


32 
649 
634 

^T 1 

634 


Galactic (GFG98) 
Galactic (FSG03) 
Galactic (T03) 


-3.037 ±0.138 

-3.06 ±0.11 
-3.141 ±0.100 


-1.021 ±0.040 
-0.989 ± 0.034 
-0.826 ±0.119 


0.209 
0.24 


28 
32 
53 


a For LMC PL relations, Assume Hlmc 


= 18.50mag. 






Table 3. The adopted I band Period-Luminosity relations. 51 




Relation 


Slope (ai) 


ZP (6/) 


ti 


N 


LMC (MF91) 
OGLE LMC (U99) 
OGLE LMC (Here) 
LMC PL >10d (TR02) 
LMC PL (Max) 


-3.06 ±0.07 
-2.962 ± 0.021 
-2.965 ±0.028 

-2.82 ±0.13 
-2.958 ±0.033 


-1.81 ±0.03 
-1.942 ±0.014 
-1.889 ±0.019 

-2.09 ±0.15 
-2.129 ± 0.023 


0.18 
0.109 
0.145 
0.12 
0.171 


32 
658 
634 
~47 
634 


Galactic (GFG98) 
Galactic (FSG03) 
Galactic (T03) 


-3.329 ±0.132 

-3.24 ±0.11 
-3.408 ±0.095 


-1.435 ±0.037 
-1.550 ±0.034 
-1.325 ±0.114 


0.194 
0.23 


27 
32 
53 



a For LMC PL relations, Assume [ilmc = 18.50ma<?. 




Fig. 1. Two examples of unusual light curves in OGLE Cepheids. Upper panels are V band light curve and lower 
panels are I band light curves. The light curves for Cepheids that have similar periods are shown with dashed lines, 
for comparison. Original data points are also indicated. 



from the OGLE web-page 4 . This consisted of 771 funda- 
mental mode Cepheids (as judged by Udalski et al. 1999b). 
Then, short period Cepheids (log(P) < 0.4) were elimi- 
nated from the sample, as in the case of Udalski et al. 



4 The URL is |http:/ /bulge.prince ton.edu/~ogle/| 



(1999a). The V and I band photometric data for the re- 
maining Cepheids were then fit with the 4*'*-order Fourier 
expansion using simulated annealing techniques (equation 
(3)) to reconstruct the light curves and obtain the mean 
magnitudes (Ngeow et al. l2003[) . The light curves for each 
Cepheid were then visually inspected, and further elimi- 
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Fig. 2. Some examples of poorly constructed V band light curves, due to bad phase coverage or clustering of data 
points in certain phases. Original data points are also indicated. 
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Fig. 3. Some examples of Cepheids that have unusually small amplitudes. The light curves for the Cepheids that have 
similar periods are shown in dashed lines, for comparison. Original data points are also indicated. 



nation of data was carried out according to the following 
criteria: 

1. In order to have a self-consistent number of Cepheids 
in both bands, Cepheids without V band data were 
removed from the sample. 

2. There are 3 Cepheids that have very unusual light 
curves, with two of them shown in Figured] The nature 
of these objects or the reason for having the unusual 
light curves is still unknown. However it is clear that 
they must be eliminated from the sample. 

3. Then, we visually examined the V band data and cor- 
responding Fourier expansions for each star to deter- 



mine which stars had acceptable Fourier expansions. 
In some cases the V band data are clustered around 
a phase point resulting in numerical bumps in the 
Fourier expansion (see Ngeow et al. p003 for more 
details). This is never a problem in the I band be- 
cause the number of points per light curve is so much 
greater (more than 120, as compare to V band which 
has ~ 12-~ 50 data points per light curves). However, 
the Cepheids with poorly reconstructed V band light 
curves were eliminated in order to have exactly the 
same Cepheids in both bands. The number of stars re- 
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Fig. 4. Two examples of well reconstructed light curves. The upper panels are V band light curves and the lower 
panels are I band light curves. Original data points are also indicated. 



jected according to this criterion is 34. Some examples 
of this case are shown in Figure |21 
4. Finally, some Cepheids appear to have a flat or low 
amplitude light curves when compared to light curves 
of stars with similar periods. Figure shows a few ex- 
amples of this case. Since we eliminated the Cepheids 
with log(P) < 0.4, these variables are unlikely to be 
the RR Lyrae-type variables. Due to the uncertainty 
of their classification, we eliminated all 43 variables 
falling into this category. We further discuss this group 
of Cepheids in Section 4.4. 

We fit the PL relations to the remaining sample (two 
examples are shown in Figure with the published ex- 
tinction values given by OGLE team, assuming Ay = 
3.24E{ B - V) and A, = 1ME{B - V) (Udalski et al. 
199QEJ), to obtain the dereddened PL relations. The ex- 
tinction map of OGLE LMC Cepheids is derived from 
observations of red clump stars along 84 lines-of-sight to- 
wards the LMC (Udalski et al. I1999bf> . Although there 
are some criticisms of the OGLE extinction map (e.g., see 
Feast EIEH Beaulieu et al. lMiTl Fouque 12003, we adopt 
the same extinction map to be consistent with the OGLE 
team. After fitting the PL relations, there are a few out- 
liers showing up in the plots (not shown) that are more 
than ~ 4cr away from the ridge lines at given periods. If 
we assume a normal distribution, these outliers can be re- 
jected (as in the case of sigma-clipping). However, these 
outliers are systematically fainter than other Cepheids at 
similar periods. For example, the V band magnitudes of 
these outliers are about l.Omag fainter than the mean 
magnitudes predicted from the ridge lines. Hence, there 
is a suspicion that these outliers could be the Type II 
Cepheids (or W Virginis), as they are generally fainter 



than the (classical) Type I Cepheids. 5 Seven outliers are 
removed according to this criterion. 

The resulting PL relations with this final sample are 
presented in Table [2] and |3J and the corresponding plots 
are given in FigureEJa). In addition, we did not fit the PL 
relations to long period Cepheids, as in Tammann et al. 
(2001), although it can be done easily. In Figure GJb), we 
also plot out the positions of the eliminated Cepheids in 
PL plots. As can be seen from the figure, most of the out- 
liers are rejected, as with sigma-clipping. However, some 
eliminated Cepheids fall along the ridge lines and would 
not be rejected by the sigma-clipping method. Our philos- 
ophy is that it is better to select bona fide Cepheids, and 
remove dubious Cepheids that can be eliminated based on 
physical grounds. Further investigation of the eliminated 
Cepheids, e.g. the low amplitude Cepheids or the three 
Cepheids with unusual light curves, are outside the scope 
of this paper. It is true that in the end there is little differ- 
ence in our V and I band PL slope relation to that initially 
published by the OGLE team and used by the KP. Never- 
theless, in the age of precision cosmology it is important 
to use appropriate data. Further our OGLE sample will be 
important when comparing with pulsation model results. 
Moreover, with mounting evidence that the PL relation in 
the LMC is " broken" , it is more appropriate to fit two PL 
relations to the LMC data. However we refit the OGLE 
data with one relation and included MF91 in our analy- 
sis to thoroughly test the effect of using our new way of 

5 The sigma-clipping algorithm will automatically remove 
these outliers, but there is no physical reasoning for doing that. 
In addition, the mean magnitudes for these stars have been cor- 
rected for extinction of Av ~ 0.5mag, unless these stars suffer 
higher extinction than the OGLE extinction map in order to 
produce the faint magnitudes. 
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Fig. 5. Left: (a) The V and I band PL relations from the final OGLE LMC sample, after correction for reddening. 
The best fit PL relations are indicated in solid lines. The dashed lines show the 2a dispersion, as expected from the 
width of instability strip. Right: (b) The locations of eliminated Cepheids in PL relations. The symbols are: filled 
circles = 3 Cepheids with unusual light curves (Figure [TJ; open circles = 34 Cepheids with poorly reconstructed light 
curves (Figure |2J; filled squares = 43 Cepheids with low amplitudes (Figure OJ) and open triangles = 7 Cepheids that 
appear dimmer (by ~ l.Omag in V band) from the ridge line at given periods. 



calculating V and I band means. This is our motivation 
for studying T02 (the LMC PL relation only for LMC 
Cepheids with periods longer than 10 days). It must also 
be pointed out that the calibrating Cepheids in the LMC 
and Galaxy have periods less than 25 days but many of 
the extra-galactic Cepheids observed by HST have peri- 
ods longer than 25 days. This is common to our work and 
that of KP and STS. 

3.2.2. The Galactic PL Relations 

Due to high extinction and less accurately known dis- 
tances to Galactic Cepheids from other independent meth- 
ods, Galactic PL relations have been considered less fa- 
vorably in the extra-galactic distance scale problem than 
LMC PL relations (e.g., Kennicutt et al. H995[) . However, 
the idea of using Galactic PL relations has been revived 
recently (Feast EES] 12303 Tammann lMiHl Fouque et al. 
2003] Thim et al. I2()()3|l because not only is the metallic- 
ity in many of the HST observed target galaxies closer to 
Galactic values, but also the current calibration of Galac- 
tic PL relations has recently been improved (e.g., see Feast 
2003 Fouque et al. 2003). In addition, because of some 
possible problems associated with the OGLE LMC PL re- 
lations, including the break at 10 days and the dominance 
of short period Cepheids (see previous subsection), and 
the possible small "depth effect" of LMC Cepheids (Groe- 
newegen 12.0001 Nikolaev et al. 2003 i. using the Galactic 
PL relations to calibrate Cepheid distances could be more 
desirable and can certainly provide an independent check 



of the LMC based scale which is independent of the LMC 
distance modulus: this is currently the largest source of 
systematic uncertainty in the extra-galactic distance scale 
(ztO.lmag). This does not imply that the systematic un- 
certainly associated with Galactic PL relations is smaller, 
but if this is true, the determination of Cepheid distances 
can be improved. 

We pick three recent Galactic PL relations which have 
both V and I band PL relations. The first Galactic PL re- 
lation is derived from 28 Cepheids using the Barnes-Evans 
surface brightness technique (Gieren et al. 119981 hereafter 
GFG98). The updated version of GFG98 PL relations has 
become available just recently. This includes 32 Galactic 
Cepheids in their sample (Fouque et al. 120031 hereafter 
FSG03). The last Galactic PL relation is adopted from 
Tammann et al. (2 0T]5|l . and includes the 28 Cepheids 
in GFG98 sample with an additional 25 Cepheids from 
Feast IJT3331 . referred as T03 in Table and We did 
not include the Galactic PL relations resulting from the 
Hipparcos calibration because these PL relations assumed 
the V band slopes are the same as in LMC PL relations, 
and hence calibrate the zero-point only (e.g., see Feast 
& Catchpole 119971 Lanoix et al. I1999bl Groenewegen & 
OudmaijerEOOni Feast EESl ■ 

3.2.3. LMC PL(Max) Relations 

Kanbur & Hendry l|1996f) . on the basis of hydrodynamical 
pulsation calculations (Simon et al. H993[) . suggested that 
PL relations at maximum light may have smaller scat- 
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Table 4. The CTE and metallicity corrections (see text for details). 
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N C eph 


[0/H] a 
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n r»7 
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HPP 1 

IN v i V lOUJ 


A<7 
~t 1 


q Qfi _i_ n on 


o no9 


finis 


07 
-u.u l 


NGC 2090 


30 


8 80 + 15 

\J . U V 7 — 1 — \ f . 1 ' } 


0.060 


-0.014 


-0.07 


NGC 2541 


29 


8.50 ±0.15 


0.0 


-0.074 


-0.07 


NGC 3031 


17 


8.75 ±0.15 


0.050 


-0.024 


0.0 C 


NGC 3198 


36 


8.60 ±0.15 


0.020 


-0.054 


-0.07 


NGC 3319 


33 


8.38 ±0.15 


-0.024 


-0.098 


-0.07 


NGC 3351 


48 


9.24 ±0.20 


0.148 


0.074 


-0.07 


NGC 3621 


59 


8.75 ±0.15 


0.050 


-0.024 


-0.07 


NGC 4321 


42 


9.13 ±0.20 


0.126 


0.052 


-0.07 


NGC 4414 


8 


9.20 ±0.15 


0.140 


0.066 


-0.07 


NGC 4535 


47 


9.20 ±0.15 


0.140 


0.066 


-0.07 


NGC 4548 


24 


9.34 ±0.15 


0.168 


0.094 


-0.07 


NGC 4725 


15 


8.92 ±0.15 


0.084 


0.010 


-0.07 


NGC 7331 


13 


8.67 ±0.15 


0.034 


-0.040 


-0.07 


STS Galaxies 


IC 4182 


27 


8.40 ± 0.20 


-0.020 


-0.094 


0.0 C 


NGC 3627 


25 


9.25 ±0.15 


0.150 


0.076 


0.05 


NGC 3982 


14 


(8.9 ±0.4) 


0.080 


0.006 


0.05 


NGC 4496A 


45 


8.77 ±0.15 


0.054 


-0.020 


0.05 


NGC 4527 


13 


(8.9 ± 0.4) 


0.080 


0.060 


0.05 


NGC 4536 


31 


8.85 ±0.15 


0.070 


-0.004 


0.05 


NGC 4639 


15 


9.00 ±0.15 


0.100 


0.026 


0.05 


NGC 5253 


5 


8.15 ±0.15 


-0.070 


-0.144 


0.0 C 


WM Galaxies 


NGC 4258 


15 


8.85 ±0.15 


0.070 


-0.004 


0.0 C 



a [O/H] = 12 + \og(0/H), adopted from Ferrarese et al. <|21TOI . except for NGC 3982 (Stetson & Gibson 200IJ and NGC 
4527 (Gibson & Stetson . 

b S z = (-0.2 ± 0.2) mag dex _1 ([0/ H] Tef - [O/H]). For LMC, [0/H] ref = 8.50dea; and for Galactic, [0/H] ref = 8.87dex. 
c Since the observations of NGC 3031, IC 4182 and NGC 5253 were done with WFPC (not WFPC2), CTE correction is not 
included. The reduction of NGC 4258 used Stetson's (1998) calibration (Newman et al. 2001), no CTE correction is needed. 



ter. However, they found that for a sample of 32 stars 
in the LMC, the PL relation at maximum light (here- 
after PL (Max)) in the LMC had comparable scatter to 
its counterpart at mean light. Motivated by this, and to 
use it primarily as a check on our mean light results, we 
computed the distances to all target galaxies studied in 
this paper using PL (Max) relations. Maximum light for 
Cepheids in the target galaxies are estimated using the 
light curve reconstruction techniques described in Ngeow 
et al. l|2TMJ> . 

We use the same OGLE LMC Cepheid sample as in 
the previous subsection to derive the PL (Max) relations. 
The estimation of maximum light (or the equivalent mini- 
mum magnitude) of a Cepheid from its light curve is more 
severely influenced by the quality of reconstructed light 
curves, as compared to its mean light counterpart. This is 
because bad phase coverage in the data will result in nu- 
merical bumps in the reconstructed light curves (Ngeow 
et al.EnnH which can be higher than the maximum light. 



Therefore, the elimination of poorly reconstructed light 
curves is more important in this aspect. After the maxi- 
mum light of each Cepheid in the sample has been esti- 
mated from the reconstructed light curves, the PL (Max) 
relations can be obtained. For this work, we assume that 
the reddening at maximum light is as same as at mean 
light even though the period color relations at maximum 
light are different to those at mean light (Simon et al. 
EOT . 

3.3. CTE and Metallicity Corrections 

The charge-transfer efficiency (CTE) for WFPC2 is a com- 
plicated issue. Simply speaking, it has been found that 
the performance of the WFPC2 depends on the exposure 
time, which lead to the "long-vs-short" exposure correc- 
tions, and other factors (see Stetson fl998l Freedman et al. 
2001 for more details). The published photometry data for 
KP galaxies were based on the Hill et al. (1998) calibra- 
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tion, but the final results presented in Freedman et al. 
I|2001[) are in the Stetson 1J1998II calibration. To convert 
the calibration from the Hill system to the Stetson sys- 
tem, a correction of Scte — — 0.07± 0.02mag is added to 
the distance modulus (Mould et al. 200D3 Freedman et al. 
2001). In contrast, the STS galaxies adopted the Holtz- 
man et al. (1995 ) calibration, which leads to a correction 
of Scte = 0.05mag (see Saha et al. 1200 lal for example). 
The calibration used in the WM galaxy is in the Stetson 
system (Newman et al. l2001|l . hence no CTE correction is 
needed. Note that no CTE correction is added to the three 
galaxies (IC 4182, NGC 3031 & NGC 5253) observed with 
WFPC in the early days of HST. The CTE corrections for 
all target galaxies are summarized in the last column of 
Table 21 and the CTE corrected distance modulus will be 
represented as Ho,cte- 

Due to metallicity differences between target galaxies 
and the LMC, metallicity corrections, as given in equa- 
tion (2), are commonly applied to /io,cT£ (Freedman et 
al. 12001(1 . The metallicity of the target galaxies, in terms 
of [O/H], is given in column 3 of Table 0J and the corre- 
sponding metallicity corrections with respect to the LMC, 
S Z (LMC), are tabulated in column 4 for the same table. 
Although the (average) metallicity of the target galaxies 
is closer to the Galactic value than the LMC (hence the 
metallicity correction is small), there is no obvious reason 
for not applying the same metallicity correction when us- 
ing a Galactic calibrating PL relation. Therefore, we apply 
the same correction, 5 Z (GAL), to the distance moduli as 
in the case of using LMC PL relations (thought see Section 
4.4). These metallicity corrections are listed in column 5 
of Table 01 for each target galaxies. 



4. Results 

We use the same Cepheids in each target galaxy that were 
used by the original KP or STS study after various selec- 
tion criteria (including the period-cut, color-cut and the 
quality of the Cepheids) had been applied. The periods 
of the Cepheids in final samples were taken from the cor- 
responding papers. The mean magnitudes were obtained 
with the Fourier techniques described in Section 3.1. Then 
the reddening-corrected distance modulus to individual 
Cepheids was calculated from equation (1) by using ei- 
ther the LMC or Galactic PL relations. The unweighted 
mean of the distance moduli to individual Cepheids was 
taken to be the final distance modulus to the target galaxy. 
The results of the distance moduli are presented in Tabled 
and|H|when using both the LMC and Galactic PL relations 
respectively. The results obtained from LMC PL (Max) re- 
lations are given in Table[71 Note that the distance moduli 
in these tables have been corrected for CTE (column 6 in 
Table 0J but not for metallicity effects, and all the errors 
in the tables are random (statistical) errors only. 



4-1. Comparisons of the Mean Magnitudes 

We compare the difference of mean magnitudes obtained 
with our method to published values. The results of this 
comparison are presented in TableEl where we give the av- 
erage difference of mean magnitudes (A/TV) in both bands 
for each target galaxy. Overall, the differences for the 16 
KP galaxies are: < AV >= 0.039 ± 0.012mag. and < 
Al >= 0.032 ±0.009ma ff .; for 8 STS galaxies: < AV > = 
0.035±0.010mag. and < Al >= -0.021±0.008maff.; and 
for all 25 target galaxies: < AV >= 0.036±0.008ma#. and 
< Al >= 0.014 ± 0.008mag. 

We also compare our distance moduli, calculated us- 
ing our Fourier techniques to estimate the V and I band 
means, with published distance moduli. The only two PL 
relations we can use in this comparison are the LMC PL 
relations from MF91 and U99, for STS and KP+WM 
galaxies, respectively. The results of these comparisons are 
presented in Table ^3 which shows the average difference 
between our distance moduli and the published values. In 
overall, the difference is small among all the target galax- 
ies, with a difference of — 0.020±0.017ma<7. (corresponding 
to ~1% change in distance). This result implies that our 
method of calculating means is a viable alternative tech- 
nique, as can be seen from Figure Ngeow et al. |2003 ) 
show that our method has advantages in some situations. 
The largest and smallest difference comes from NGC 4639 
and NGC 4414, with a difference of -0.240 ± 0.245mag. 
and 0.004 ± 0.066ma<?., respectively. 

While STS use phase weighted averages to compute 
their means, the KP use light curve template techniques 
(Stetson [1996 ) . The difference between this and the 
method adopted here is described in detail in Ngeow et 
al. ( 2003 ). Both methods fit equations like equation (3) to 
the observed data. In the case of Stetson l|1996fl . Ai is ob- 
tained from the data while Ak , k > 1 is obtained from the 
templates showing how the ratio of Ak/Ai,k > 1 varies 
with period. In our case all four parameters Ak ,k>l 
are obtained from the data with a simulated annealing 
approach. We do this, firstly, because we want to recon- 
struct the light curve to model it and, secondly, to capture 
any light curve shape changes due to metallicity. The top 
two panels of Figured show the second Fourier amplitude 
coefficient (A 2 ) calculated with this method for very well 
sampled Galactic (left panel) and LMC (right panel)data. 
Since this is very well sampled data, it is safe to assume 
that the progression of A 2 with period is a good represen- 
tation of reality. The bottom panel shows the same Fourier 
amplitude calculated using the Stetson (1996} method. It 
can clearly be seen that at all periods the range of A 2 at 
given period is larger than would be predicted by the Stet- 
son l)1996|l technique. Furthermore, from the well sampled 
data, we randomly pick 12 points (so that they could be 
clustered or not as the case may be), add Gaussian errors 
to these points and then calculate a simulated annealing 



12 Kanbur et al.: Cepheid distance from LMC and Galactic PL relations 

Table 5. The Cepheid distance (with CTE correction) to nearby galaxies with different LMC PL relations. 
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NGC 2541 
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071 
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348 
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071 
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27 
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27 


791 
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27 
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093 


NGC 3198 
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30 
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30 


721 
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30 
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060 


NGC 3319 


30 


593 


± 
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30 


511 


± 





092 


30 


501 


± 





092 


30 


528 ±0 
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NGC 3351 


29 


979 


± 





090 


29 


914 


± 





090 


29 


902 


± 





090 


29 


927 ±0 


090 


NGC 3621 


29 
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± 





059 


29 
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± 
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29 
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± 
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NGC 4321 
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± 





067 


30 


769 


± 





067 


30 


764 


± 





067 


30 


798 ±0 


067 


NGC 4414 


31 


236 


± 





052 


31 


104 


± 





043 


31 


099 


± 





044 


31 


133 ±0 


045 


NGC 4535 


31 


001 


± 





054 


30 


879 


± 





053 


30 


874 


± 





053 


30 


906 ±0 


053 


NGC 4548 


30 


867 


± 





076 


30 


783 


± 





078 


30 


773 


± 





078 


30 


800 ±0 


078 


NGC 4725 


30 


522 


± 





066 


30 


399 


± 





067 


30 


393 


± 





067 


30 


426 ±0 


067 


NGC 7331 


30 


891 


± 





083 


30 


800 


± 





078 


30 


791 


± 





079 


30 


819 ±0 


079 


IC 4182 


28 


365 


± 





128 


28 


334 


± 





130 


28 


318 


± 





130 


28. 


339 ± 0. 


129 b 


NGC 3627 


30 


238 


± 





083 


30 


134 


± 





084 


30 


127 


± 





084 


30 


157 ±0 


083 


NGC 3982 


31 


714 


± 





136 


31 


599 


± 





134 


31 


593 


± 





134 


31 


624 ±0 


135 


NGC 4496A 


30 


959 


± 





043 


30 


839 


± 





043 


30 


833 


± 





043 


30 


865 ±0 


043 


NGC 4527 


30 


755 


± 





129 


30 


635 


± 





127 


30 


629 


± 





128 


30 


661 ±0 


128 


NGC 4536 


30 


961 


± 





102 


30 


830 


± 





102 


30 


826 


± 





102 


30 


859 ±0 


102 


NGC 4639 


31 


790 


± 





108 


31 


664 


± 





105 


31 


659 


± 





106 


31 


692 ±0 


106 


NGC 5253 


27 


881 


± 





244 


27 


873 


± 





246 


27 


854 


± 





246 


27. 


872 ± 0. 


245 c 


NGC 4258 


29 


435 


± 





058 


29 


384 


± 





056 


29 


370 


± 





056 


29 


393 ±0 


057 



a MF91 = Madore & Freedman (I199H PL relations; U99 = Udalski et al. 11999at PL relations; Here = PL relations derived 
in Section 3.2.1; TR02= Tammann & Reindl PS) PL relations for log(P) > 1.0. See Section 3.2.1 for details. 
b There are 12 Cepheids in this galaxy with period less than 10 days. If we use the PL<iod relations, as given by Tammann 
& Reindl (2002J, then fi — 28.340 ± 0.130, which is identical to the value with PL>iod relations. 

c There are 2 Cepheids in this galaxy with period less than 10 days. If we use the PL<iod relations, as given by Tammann 
& Reindl ( 2002), then fi = 27.876 ± 0.246, which is very close to the value with PL>iod relations. 



fit, we recover trends in the top panel of Figured This is 
clearly displayed in figure 18 of Ngeow et al. ( 2003 ). 



4-2. Distance Moduli from Various PL Relations 

Since the distance moduli listed in Tabled are calculated 
using different LMC PL relations, they can be compared. 
A quick comparison in Table d indicates that the MF91 
PL relation always give the longer distance moduli, and 
the shortest distance moduli are from the OGLE LMC PL 
relations derived in Section 3.2.1. The average difference 
between the distance moduli derived from MF91 and the 
OGLE LMC (Here), for all target galaxies in Tabled is 
about 0.105 ± 0.006mag. Comparing this result with the 
previous case (i.e. Section 4.1) of using the same PL rela- 
tions but different means magnitudes, the distance modu- 
lus is more sensitive to different calibrating PL relations. 
Also, the distance moduli from three OGLE LMC PL rela- 
tions are consistent with each other within the statistical 
errors. 



A similar situation exists when using Galactic PL rela- 
tions, as shown in Tabled The distance moduli from both 
GFG98 and FSG03 Galactic PL relations are consistent 
with each other (with a difference of ~ 0.024 ± 0.007ma<7., 
as there are 26 common Cepheids in both samples), while 
FSG03 PL relations produce a shorter distance. However, 
the distance moduli from T03 PL relations are systemat- 
ically further than the other two, although all three PL 
relations share some common Galactic Cepheids. Never- 
theless, the average difference between the distance moduli 
from T03 and FSG03 is ~ 0.097 ± 0.008ma 5 ., comparable 
to the difference seen in LMC PL relations. 

Finally, we compare the distance moduli (after CTE 
corrections) from LMC PL(Max) relations in Table d to 
their mean light counterparts, i.e., column 4 of Tabled 
The distance moduli from PL (Mean) and PL (Max) rela- 
tions are consistent with each other, although PL(Max) 
relations generally give a shorter distance. This serves 
as an important check on our calibration of the extra- 
galactic distance scale. The closer distance from PL (Max) 
relations is expected (with the difference of ~ 0.015 ± 
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Table 6. The Cepheid distance (with CTE correction) to nearby galaxies with different Galactic PL relations. 



13 



Galaxy 



NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 



925 

1326A 

1365 

2090 

2541 

3031 

3198 

3319 

3351 

3621 

4321 

4414 

4535 

4548 

4725 

7331 



IC 4182 
NGC 3627 
NGC 3982 
NGC 4496A 
NGC 4527 
NGC 4536 
NGC 4639 
NGC 5253 



NGC 4258 



GFG98 a 



FSG03 a 



T03 a 



29.765 
31.093 
31.336 
30.415 
30.464 
27.867 
30.842 
30.569 
29.936 
29.249 
30.929 
31.265 
31.020 
30.846 
30.541 
30.875 



±0.061 
± 0.090 
± 0.053 
± 0.038 
± 0.075 
± 0.098 
±0.061 
± 0.096 
± 0.090 
± 0.060 
±0.067 
± 0.063 
± 0.055 
± 0.078 
±0.067 
±0.091 



29.768 ± 0.060 
31.070 ± 0.089 
31.283 ± 0.052 
30.392 ± 0.036 
30.430 ± 0.072 
27.854 ± 0.095 
30.805 ± 0.060 
30.560 ± 0.094 
29.946 ± 0.090 
29.226 ± 0.059 
30.869 ± 0.067 
31.204 ±0.052 
30.969 ± 0.054 
30.834 ± 0.078 
30.490 ± 0.066 
30.858 ± 0.083 



29.833 ±0.061 
31.166 ±0.091 
31.414 ±0.053 
30.488 ± 0.039 
30.539 ± 0.075 
27.939 ± 0.098 
30.917 ±0.061 
30.639 ± 0.096 
30.004 ± 0.090 
29.322 ±0.061 
31.008 ±0.067 
31.344 ± 0.064 
31.097 ±0.055 
30.917 ± 0.078 
30.618 ± 0.067 
30.947 ± 0.092 



28.285 
30.237 
31.726 
30.974 
30.771 
30.989 
31.812 
27.776 



±0.129 
± 0.083 
±0.139 
± 0.043 
±0.132 
±0.102 
±0.112 
± 0.243 



28.330 
30.206 
31.682 
30.926 
30.723 
30.929 
31.758 
27.846 



±0.128 
± 0.083 
±0.136 
± 0.043 
±0.129 
±0.102 
±0.109 
± 0.243 



28.346 ± 
30.311 ± 
31.802 ± 
31.051 ± 
30.848 ± 
31.068 ± 
31.890 ± 
27.834 ± 



0.129 
0.083 
0.139 
0.044 
0.132 
0.102 
0.113 
0.242 



29.377 ± 0.062 29.401 ± 0.058 29.442 ± 0.062 



a GFG98 = Gieren et al. flTOt PL relations; FSG03 = Fouque et al. (I2U031 PL relations; T03 = Tammann et al. (I2U031 
PL relations. See Section 3.2.2 for details. 



O.OOAmag., compared to the mean light counterparts), be- 
cause Cepheids will appear closer at maximum light. The 
errors associated with maximum light are generally com- 
parable to those at mean light. This will be studied in 
future work. 

4-3. Comparison to Published Results 

In order to compare our results with published values, 
we need to select the representative distance moduli from 
both LMC or Galactic PL relations. We pick the distance 
moduli from the derived OGLE LMC PL relations in Sec- 
tion 3.2.1 (i.e. column 4 of Table and FSG03 PL rela- 
tions to represent the LMC and the Galactic results, re- 
spectively. The CTE-corrected distance moduli, Ho,Cte, 
are listed in Table [SJ However, the metallicity corrections 
are not required in these comparisons because they are 
the same for published results and this work. The effect 
of metallicity corrections are discussed in Section 4.4. 

The results of the comparisons are plotted in Figure 
[S] and the average difference between our distance mod- 
uli and the published values are presented in Table ITT1 for 
both the LMC and Galactic PL relations. The overall com- 
parisons indicate that our results are consistent with the 
published results. When using the LMC PL relation, our 



distance moduli agree well with the KP results, but show 
a large discrepancy for the STS results (— 0.178mag.). A 
careful comparison of Table|Slwith published results shows 
that this difference arises from two sources: ~ O.lOmag. 
is due to the use of MF91 LMC PL relations by STS and 
~ 0.07mag. is due to differences in the V and I band 
means when calculating the distance moduli (see Section 
4.1 and Table E3>- However, when using the Galactic PL 
relations, our results are different by ~ O.lmag. compared 
to both KP and STS results. Clearly, the use of Galactic 
PL relations increases the distance moduli. Surprisingly, 
the distance modulus of NGC 4258 calculated from LMC 
PL relations is closer to the water maser distance although 
the metallicity of this galaxy is nearly identical to Solar. 

Gibson et al. (2000, hereafter GOO) used KP techniques 
to reanalyze the STS galaxies. Since GOO did the photom- 
etry from scratch, their list of Cepheids in each target 
galaxy is in general different from those found by STS. 
Moreover, they used the Madore & Freedman (Q~991) cal- 
ibrating PL relations (MF91). Gibson & Stetson lj2T)I7T)l 
and Stetson & Gibson (200 1 i used the same photometry as 
GOO but used the Udalski et al. l)1999a|l LMC PL relations 
(U99) to re-calibrate the distance to all the STS galaxies. 
They found shorter distance than those published by STS. 
They suggested that the discrepancy between their work 
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Fig. 6. Upper panel: Comparison 
of the distance moduli (after cor- 
rected for CTE) obtained by us- 
ing the Fourier techniques to es- 
timate the mean magnitudes with 
the published values. The line rep- 
resents the case of iip u b,CTE = 
H*Here,CTE- Lower panel: The plot 
of the difference, A = HHere.CTE — 
I^Pub,CTE, vs. the distance moduli. 
The dashed line is the average dif- 
ference of —0.02. The symbols are: 
crosses = 16 KP galaxies; open cir- 
cles = 8 STS galaxies; and filled cir- 
cle = WM galaxy. 



and that of STS is due to the use of U99 LMC PL re- 
lations, their photometric analysis of the raw data and 
the STS analysis to account for possible correlated mea- 
surement errors in V and I. While these are important 
issues, our approach is to use exactly the same photome- 
try as used by both KP and STS groups. Our methods are 
closer to KP, though we have a different way of estimating 
the V and I band means. Our results show that, even with 
exactly the same Cepheids as used by STS, our distance 
moduli are significantly shorter than STS but significantly 
longer than those published in Gibson & Stetson (|2001). 
It would be interesting to perform our analysis on the GOO 
photometry but this is not available as advertised on the 
KP web-sitc. 



4-4- Comparison of Metallicity- Corrected Distance 
Moduli 

As well as comparing our results with the published val- 
ues, we discuss the distance moduli obtained when LMC 
and Galactic PL relations are used. First, comparing 
the distance moduli that are uncorrected for metallicity, 
Ho,cte{LMC) and /jo,cte(GAL) in Table show that 
the distance moduli from LMC PL relations are always 
shorter than their Galactic counterparts (see FigureEfa)). 
Summaries of this comparison are given in Table IT2"1 This 



shows that the LMC PL relations will give a smaller dis- 
tance modulus by ~ 0.07 ± O.Olmag on average (or about 
3.5% in distance) compared to the distance modulus ob- 
tained from Galactic PL relations. The negligible random 
errors suggest that the difference of 0.07 'mag. might be 
significant, although it is small. However, when the metal- 
licity corrected distance moduli from the two sets of PL 
relations (fj, z (LMC) and fi z (GAL) in Table© are com- 
pared, the difference between them falls close to zero. The 
results are listed in Table IT21 and shown in Figure EJb). 

Another approach is to use the Galaxy and LMC 
(P> lOd) as calibrating PL relations for metal rich and 
metal poor galaxies respectively. Thus if the seven SN 
calibrators of STS were forced on an LMC PL relation 
(U99), one would obtain an average distance decrease of 
0.17 'mag. However, if we force the five galaxies in their 
sample, which have on average the same metallicity as 
the Galaxy, onto the steep relation given in T03, the pub- 
lished STS distances are recovered for these galaxies to 
within O.Olmag. Applying the same procedure to the 10 
KP galaxies which are metal rich increases their distance 
moduli by about 0.16mag. on average, whilst using TR02 
for the remaining metal poor galaxies increases their dis- 
tance moduli by 0.02mag. 

It has been pointed out by an anonymous referee 
that it may not be appropriate to eliminate those 
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0.5 



iog(P) 



log(P) 



ioa(P) 



loS(P) 



Fig. 7. Comparison of Ai from Fourier techniques and template methods. The upper panels show the values of Ai 
from direct Fourier fit with the method described in Ngeow et al. (2002J), an( i the lower panels show the values of 
A2 by using the template light curves given by Stetson (1996). Since the template light curves are undefined for 
log(p) < 0.85, Cepheids with period shorter and longer than log(p) = 0.85 are represents as open circles and filled 
triangles, respectively. Left:(a) The "calibrating set" Cepheids. Rightih) The OGLE LMC Cepheids. 



Table 7. The Cepheid distance to nearby galaxies with 
PL (Max) relations. 



Galaxy 


OGLE LMC PL(Max) a 


NGC 925 


29.699 


± 


0.059 


NGC 1326A 


30.994 


± 


0.089 


NGC 1365 


31.165 


± 


0.051 


NGC 2090 


30.309 


± 


0.035 


NGC 2541 


30.343 


± 


0.071 


NGC 3031 


27.793 


± 


0.094 


NGC 3198 


30.712 


± 


0.062 


NGC 3319 


30.479 


± 


0.092 


NGC 3351 


29.871 


± 


0.087 


NGC 3621 


29.147 


± 


0.059 


NGC 4321 


30.740 


± 


0.066 


NGC 4414 


31.099 


± 


0.042 


NGC 4535 


30.856 


± 


0.053 


NGC 4548 


30.750 


± 


0.081 


NGC 4725 


30.388 


± 


0.066 


NGC 7331 


30.787 


± 


0.082 


IC 4182 


28.253 


± 


0.130 


NGC 3627 


30.116 


± 


0.087 


NGC 3982 


31.592 


± 


0.143 


NGC 4496A 


30.840 


± 


0.044 


NGC 4527 


30.616 


± 


0.128 


NGC 4536 


30.843 


± 


0.102 


NGC 4639 


31.659 


± 


0.112 


NGC 5253 


27.773 


± 


0.226 


NGC 4258 


29.352 ± 0.057 



a Corrected for the CTE effects of Scte (taken from last 
columns of Table 2J. 



stars with a low pulsational amplitudes shown in Fig- 
ure |21 as this may bias the LMC PL relation. Such 
low amplitude Cepheids do exist and are seen in ex- 
ternal galaxies. If we include this group of stars in 
deriving the LMC PL relations, the new PL relations 
are: M v = -2.713(±0.044) log(P) - 1.418(±0.030) and 
Mi = -2.943(±0.029)log(P) - 1.899(±0.020). The dis- 
tance moduli derived from these new LMC PL relations is 
~ O.Olmag. shorter when compared to the distance mod- 
uli derived from LMC PL relations given in Section 3.2.1. 
When comparing the distance moduli from these new 
LMC PL relations to those derived using Galactic PL rela- 
tions, the average difference for the 25 target galaxies con- 
sidered here is -0.052±0.006ma#. and -0.082±0.006ma5. 
with and without the metallicity corrections respectively. 

These low amplitude stars are not first overtone pul- 
sators since their light curve Fourier parameters fall in 
the region occupied by fundamentals. They lie in the in- 
stability strip and do not have unusual colors for their 
period. Thus their luminosities, masses and temperatures 
are similar to other Cepheids of similar period. So why 
do they have such a low amplitude? It could be that they 
have a slightly different composition or are just entering 
or leaving the fundamental mode instability strip (Buch- 
ler & Kollath 2002). In this case it is our contention that 
they should be excluded from the sample since they are 
Cepheids undergoing a transition. 

The use of Galactic PL relations to calibrate the 
Cepheid distance scale has been tried before by Paturel 
et al. (2002a) and Paturel et al. (2002b). The first paper, 
Paturel et al. (j2D02a), used the GFG98 sample and the 
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method of "sosie" (Paturel I1984[l to determine Cepheid 
distances without assuming the PL relations. The sec- 
ond paper, Paturel et al. (2002b), applied the Hippar- 
cos calibrated PL relations (Feast & Catchpole 1997) to 
find Cepheid distances. The mean offsets of their results 
and the published KP results are 0.161 ± 0.029mag. and 
0.027 ± O.QWmag. for their first and second papers, re- 
spectively. The latter result is comparable to our work. 
However, the extra-galactic Cepheids used in both papers 
are slightly different than either the KP or STS groups, in 
contrast with our study in this paper, because these au- 
thors applied a different method to deal with incomplete- 
ness bias. Also, these authors did not apply any metallicity 
corrections in their papers. 

5. Conclusion and Discussion 

In this paper we recalibrate the Cepheid distance to about 
two dozen nearby galaxies, including the KP and STS 
galaxies, and a water maser galaxy. We use much of 
the same methodology as the KP team: the same LMC 
distance, value of R and metallicity correction, and the 
same (number of) Cepheids along with the published pho- 
tometries and periods. However our approach is differ- 
ent from the KP team in two aspects: (a) we estimate 
the mean magnitudes of sparsely sampled HST data from 
Fourier techniques (Ngeow et al. 2003) by reconstructing 
the Cepheid light curves; and (b) we use different sets of 
PL relations, including the new OGLE LMC PL relations, 
the Galactic PL relations and the LMC PL (Max) rela- 
tions, in distance determination. Overall, our results are 
consistent with each other and KP. We find significantly 
shorter distances to the STS galaxies. The use of the new 
Fourier techniques to obtain V and I band means does not 
produce significant deviations from the existing methods. 
However the derived distance modulus is more sensitive 
to the calibrating PL relation that is adopted. The most 
striking result from this study is that the distance moduli 
derived from using the LMC and Galactic PL relations are 
indistinguishable after the metallicity corrections are ap- 
plied. This provides strong support for the size and quan- 
tity of the metallicity dependence of the Cepheid PL rela- 
tions and the non-universality of PL relations (Tammann 
et al. l2UU3|> . 

Recent work has shown that the PL relation in the 
Galaxy and LMC have different slopes and moreover that 
the PL relation in the LMC is "broken" at a period of 10 
days. Metallicity affects the mean brightness of a Cepheid 
and up to now, this has been the standard motivation 
for deriving and applying simple additive corrections to 
the PL relation to account for metallicity differences be- 
tween calibrating and target galaxy. However it may be 
that metallicity also affects the slope of the PL relation. 
If this is so, then the use of a simple additive correction 
as given in equation (2) is not really appropriate because 
the PL relation may have a different slope in the target 



and calibrating galaxy. However, the right panels of figure 
8 show clearly that it works in the sense that the distance 
moduli differences obtained when using the Galactic and 
LMC as calibrating PL relations are driven close to zero. 
We try to investigate why this occurs in what follows. 

If we want to compare the distance modulus obtained 
using LMC and Galactic calibrating PL relations (or any 
pair of PL relations) when using exactly the same (number 
of) Cepheids in the target galaxy, including the same peri- 
ods and mean magnitudes, then the difference in distance 
modulus can be expressed as: 

A Mo = W - l)Aoy ~ RAaj] S y ( 5 ) 
+(R-l)Ab v - RAb h 

where Aa/y n and Abiy,i) & r e the differences in slopes 
and zero-points for the two PL relations, respectively. 
Then, for these two PL relations, the change in the dis- 
tance modulus is a simple linear function of the target 
galaxy period distribution (< log(P) >= J^log(P)i/iV, 
hereafter mean period) , under the assumption of constant 
R (i.e. the universality of Galactic extinction law, see ob- 
servational verification for this assumption by Macri et 
al. E001). The error in A/xq is the quadrature sum of 
the from both PL relations for the target galaxy. If 
metallicity corrections of 5 Z (from equation (2)) are ap- 
plied to equation (5) , the metallicity of the target galaxies 
cancels out and leaves the difference between the LMC 
and Galactic metallicity. Adopting [0/H]lmc = 8.50dea; 
and [0/H]gal — 8.87dex, along with 7 = —0.2 ± 
0.2mag. dex~ x , the difference in metallicity-corrected dis- 
tance modulus becomes: 

Afj, z = Afi - 0.2(8.50- 8.87) (6) 
= A^ + 0.074 (±0.074) 

It is worth pointing out that equations (5) and (6), 
though straightforward to derive, have not been presented 
in the literature before to the best of our knowledge. 
Equation (6) is of course dependent on the metallicity 
correction law adopted and the slopes of the calibrat- 
ing PL relations. We can re-write equation (5) as A/io = 
c < log(P) > +d, where A^o = Ho(LMC) - [i (GAL), 
c = (R-l)Aa v -RAai and d = (P- l)A6 y - PA6/. For 
the four LMC PL relations (excluding the PL (Max) re- 
lations) and the three Galactic PL relations considered 
in this paper, there are a total of 12 combinations of 
(LMC, GAL) PL relations. The coefficients of c and d 
for each combination are listed in column 2 & 3 of Table 
1131 In the same table, we also list out the coefficients of 
d after applying the metallicity corrections (i.e. equation 
(6)), i.e. A[i z =c< log(P) > +d z where d z = d + 0.074. 

To see what values of < log(P) > would produce 
identical distance moduli from using either the LMC or 
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Table 8. Comparisons of the distance moduli to target galaxies after CTE and metallicity corrections. 



Galaxy 


fJ-0,CTE 


(LMC) 


fi z (LMC) 


Ho,cte(GAL) 


fj, z (GAL) 


^cte{PUBY 


H.{PUBY 


NGC 925 


29.718 


± 


0.060 


29.728 


29.768 


± 


0.060 


29.704 


29.80 ±0.04 


29 


81 


NGC 1326A 


30.998 


± 


0.088 


30.998 


31.070 


± 


0.089 


30.996 


31.04 ±0.10 


31 


04 


NGC 1365 


31.186 


± 


0.051 


31.287 


31.283 


± 


0.052 


31.301 


31.18 ±0.05 


31 


27 


NGC 2090 


30.319 


± 


0.036 


30.379 


30.392 


± 


0.036 


30.378 


30.29 ±0.04 


30 


35 


NGC 2541 


30.348 


± 


0.071 


30.348 


30.430 


± 


0.072 


30.356 


30.25 ±0.05 


30 


25 


NGC 3031 


27.791 


± 


0.093 


27.841 


27.854 


± 


0.095 


27.830 


27.75 ±0.08 


27 


82 


NGC 3198 


30.721 


± 


0.060 


30.741 


30.805 


± 


0.060 


30.751 


30.68 ±0.08 


30 


70 


NGC 3319 


30.501 


± 


0.092 


30.477 


30.560 


± 


0.094 


30.462 


30.64 ± 0.09 


30 


62 


NGC 3351 


29.902 


± 


0.090 


30.050 


29.946 


± 


0.090 


30.020 


29.85 ± 0.09 


30 


00 


NGC 3621 


29.155 


± 


0.059 


29.205 


29.226 


± 


0.059 


29.202 


29.08 ±0.06 


29 


11 


NGC 4321 


30.764 


± 


0.067 


30.890 


30.869 


± 


0.067 


30.921 


30.78 ±0.07 


30 


91 


NGC 4414 


31.099 


± 


0.044 


31.239 


31.204 


± 


0.052 


31.270 


31.10 ±0.05 


31 


24 


NGC 4535 


30.874 


± 


0.053 


31.014 


30.969 


± 


0.054 


31.035 


30.85 ±0.05 


30 


99 


NGC 4548 


30.773 


± 


0.078 


30.941 


30.834 


± 


0.078 


30.928 


30.88 ±0.05 


31 


05 


NGC 4725 


30.393 


± 


0.067 


30.477 


30.490 


± 


0.066 


30.500 


30.38 ±0.06 


30 


46 


NGC 7331 


30.791 


± 


0.079 


30.825 


30.858 


± 


0.083 


30.818 


30.81 ±0.09 


30 


84 


IC 4182 


28.318 


± 


0.130 


28.298 


28.330 


± 


0.128 


28.236 


28.36 ± 0.09 b 






NGC 3627 


30.127 


± 


0.084 


30.277 


30.206 


± 


0.083 


30.282 


30.22 ±0.12 






NGC 3982 


31.593 


± 


0.134 


31.673 


31.682 


± 


0.136 


31.688 


31.72 ±0.14 






NGC 4496A 


30.833 


± 


0.043 


30.887 


30.926 


± 


0.043 


30.906 


31.03 ±0.14 






NGC 4527 


30.629 


± 


0.128 


30.709 


30.723 


± 


0.129 


30.783 


30.72 ±0.12 






NGC 4536 


30.826 


± 


0.102 


30.896 


30.929 


± 


0.102 


30.925 


31.10 ±0.13 






NGC 4639 


31.659 


± 


0.106 


31.759 


31.758 


± 


0.243 


31.784 


32.03 ±0.22 






NGC 5253 


27.854 


± 


0.246 


27.784 


27.846 


± 


0.243 


27.702 


28.08 ± 0.20 c 






NGC 4258 


29.370 


± 


0.056 


29.440 


29.401 


± 


0.058 


29.397 


29.40 ± 0.09 


29 


47 



a Published distance moduli for KP galaxies are taken from Freedman et al. I|200ip . Table 4. The distance moduli for STS 
galaxies and WM galaxy are taken from series of STS papers and Newman et al. ( 2001), respectively. 
b Saha et al. (|T5943 assumed Av = Aj = 0, hence fi — (fj,v + /Ur)/2. 

c The paper (Saha et al. 1995 I did not list out the final n, hence this value is calculated via equation (2) with [iv and hi 
given in the paper. 




28 29 30 31 32 28 29 30 31 



Fig. 8. Comparisons with published results. The symbols are same as in Figure|S| Left: (a) Comparison of our results 
from the new OGLE LMC PL relations. The dashed line is the average difference of —0.06. Right:(b) Comparison of 
our results from the Galactic PL relations. The dashed line is the average difference of 0.01. 
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Fig. 9. Comparisons of distance moduli from LMC and Galactic PL relations. The symbols are same as in Figure 
El Left:(a) Uncorrected for metallicity. The dashed line is the average difference of —0.07. Right:(b) Corrected for 
metallicity. The dashed line is the average difference of —0.00. 



Table 9. Comparisons of the mean magnitudes with pub- 
lished results. 



Galaxy 


N 


AV/N a 


AI/N a 


NGC 925 


72 


-0.046 ± 0.004 


-0.006 ± 0.010 


NGC 1326A 


15 


0.040 ± 0.009 


0.007 ± 0.006 


NGC 1365 


47 


-0.077 ± 0.006 


-0.029 ± 0.014 


NGC 2090 


30 


0.079 ± 0.004 


0.064 ± 0.007 


NGC 2541 


29 


0.113 ±0.015 


0.086 ± 0.020 


NGC 3031 


16 


0.065 ± 0.022 


0.039 ± 0.014 


NGC 3198 


36 


0.001 ±0.059 


0.020 ± 0.063 


NGC 3319 


33 


0.046 ± 0.006 


-0.001 ±0.012 


NGC 3351 


48 


0.035 ± 0.029 


0.009 ± 0.012 


NGC 3621 


59 


0.062 ± 0.029 


0.072 ± 0.008 


NGC 4321 


42 


0.018 ±0.009 


-0.005 ±0.013 


NGC 4414 


8 


0.027 ± 0.007 


0.077 ± 0.081 


NGC 4535 


47 


0.066 ± 0.008 


0.051 ± 0.009 


NGC 4548 


24 


0.087 ±0.021 


0.023 ±0.012 


NGC 4725 


15 


0.024 ±0.011 


0.023 ± 0.008 


NGC 7331 


13 


0.086 ±0.010 


0.077 ± 0.031 


IC 4182 


26 


0.050 ±0.013 


-0.011 ± 0.020 


NGC 3627 


25 


0.045 ± 0.008 


-0.018 ±0.026 


NGC 3982 


14 


-0.018 ±0.014 


-0.030 ± 0.014 


NGC 4496A 


45 


0.043 ± 0.004 


-0.016 ±0.019 


NGC 4527 


13 


-0.013 ±0.016 


-0.010 ±0.020 


NGC 4536 


31 


0.057 ± 0.005 


-0.005 ±0.015 


NGC 4639 


15 


0.041 ± 0.009 


-0.007 ± 0.017 


NGC 5253 


5 


0.026 ± 0.007 


-0.071 ± 0.045 


NGC 4258 


15 


0.049 ±0.018 


0.007 ± 0.012 



Table 10. Comparisons of the distance moduli with dif- 
ferent mean magnitudes 51 . 



Case 


A(mag.) b 


16 KP + WM galaxies 


0.005 ±0.015 


8 STS galaxies 


-0.075 ± 0.037 


All 25 galaxies 


-0.020 ±0.017 



a By using the same PL relations but different mean mag- 
nitudes. See Section 4.1 for details. 

b Mean Difference, A =< no,cTE{Here)—no,cTE{Pub.) >, 
the errors are the standard deviations of the means. 



Table 11. Comparisons with the published results. 



Case 






A(mag.) a 






Usin 


g LMC PL relations 


16 KP 
8 STS 
All 25 


galaxies 
galaxies 
galaxies 




-0.002 ±0.017 
-0.178 ±0.039 b 
-0.059 ±0.023 






Using 


Galactic PL relations 


16 KP 
8 STS 
All 25 


galaxies 
galaxies 
galaxies 




0.075 ±0.018 
-0.108 ±0.038 
0.014 ±0.024 



b AV = V He 



a Mean Difference, A =< fio,CTE(Here) — fj,o,cTE{Pub.) >, 
the errors are the standard deviations of the means. 
b This difference is mainly due to the different LMC PL 
relations used. See text (Section 4.3) for details. 



Same for AI. 



columns 5 & 6 in Table ^] for the case of AfiQ = 
the Galactic PL relations, we solve for < log(P) > such and A[i z = 0, respectively. From Table El the mean pe- 
that A/i = or Afj, z — 0. The solutions are listed in riod for two distance moduli to be identical is around 
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Table 12. Comparisons of metallicity-corrected distance 
moduli. 



Case 


A(mag.) a 


Without metallicity correction 


16 KP galaxies 


-0.077 ±0.005 


8 STS galaxies 


-0.070 ±0.015 


All 25 galaxies 


-0.073 ± 0.006 


With metallicity correction 


16 KP galaxies 


-0.002 ± 0.005 


8 STS galaxies 


-0.003 ±0.018 


All 25 galaxies 


-0.001 ± 0.007 



a Mean Difference, A =< fi(LMC) - /i(GAL) >, the errors 
are the standard deviations of the means. 



< log(P) > ~ 1.1 - 1.2, and < log(P) > z ~ 1.3 - 1.4, 
except for the (MF91.GFG98) and (MF91,FSG03) pairs. 
In Table El we list out the observed mean periods for 
the target galaxies. This shows that most of the extra- 
galactic Cepheids in our target galaxies have an observed 
mean period of ~ 1.4 (the mean period for 16 KP galaxies 
is 1.418; the mean period for 8 STS galaxies is 1.401; and 
the mean period for all 25 galaxies is 1.342). Therefore, 
without applying the metallicity correction, the Galactic 
PL relations will be expected to produce longer distances 
than the LMC PL relations in most of the (LMC, GAL) 
pair of PL relations. For example, the change of the dis- 
tance modulus is ~ O.lOmag. when comparing the LMC 
PL relations (either U99 or the PL relations derived in 
Section 3.2.1) to the FSG03 Galacic PL relations in in- 
dividual target galaxies. However, the LMC and Galactic 
PL relations will produce almost identical distance mod- 
uli after a metallicity correction to within ~ 0.03ma<?., for 
most of the target galaxies. The only exception is NGC 
5253, because the mean period for this galaxy is 1.029, 
which is much smaller than the required < log(P) > of 
- 1.4. 

The T03 Galactic PL relation has a different slope to 
the GFG03 and FSG03 Galactic PL relations. Thus the 
median period in Table IT3*I required for a simple additive 
correction to be sufficient is approximately 1.2. This is 
slightly outside the range of mean period required to pro- 
duce similar distance moduli after a simple metallicity 
correction when both Galactic and LMC calibrating PL 
relations are used. 

If the result of this paper, i.e. the near identical dis- 
tance moduli from LMC and Galactic PL relations af- 
ter metallicity corrections, is true, then this result can 
be used to constrain the sign for the coefficient of metal- 
licity correction, the 7 in equation (2). The value of 
7 = 0.2±0.2TOa<7 dex~ 1 used in this study is adopted from 
Freedman et al. (2001), which is roughly the midrange 
value from several empirical studies (see Section 3). Since 
the difference of distance moduli between LMC and Galac- 
tic PL relations is about —0.07mag. without the metallic- 



ity correction fTable lT2*|l . a +0.07mag. metallicity correc- 
tion is required to bring the two distance moduli to be 
identical. The correction of 0.074ma<?. in equation (6) is 
almost identical to this requirement, hence the value of 7 
should be around —0.2mag dex~ l and constrains the sign 
to be negative. In addition, if the Cepheid PL relations 
do indeed depend on metallicity, the result of this paper 
suggests that a simple additive metallicity correction as 
in equation (2) is a good approximation to model the full 
complexity of the metallicity dependence of the Cepheid 
PL relation, provided the mean period of Cepheids in 
the target galaxy are in the appropriate range for the 
slopes and extinction laws adopted for the calibrating PL 
relations 6 . Some researchers suggest using the LMC PL 
relations and Galactic PL relations for metal-poor and 
metal-rich galaxies, respectively, forgetting about metal- 
licity corrections until an solid understanding of this topic 
is obtained. This naturally begs the question of which cal- 
ibrating PL relation to use if the metallicity of the target 
galaxy is in between the LMC and Galaxy. However our 
result does not, at the moment, provide evidence support- 
ing one Galactic PL relation over another. 

In summary, the above discussion suggests that: (a) 
metallicity corrections are necessary when using Cepheid 
PL relations to find distance moduli; (b) as a consequence, 
the PL relations do depend on metallicity; (c) hence, 
Cepheids in the LMC and Galaxy obey different PL rela- 
tions; (d) the sign for the coefficient of metallicity correc- 
tion (7) has to be negative; and (e) both LMC and Galac- 
tic PL relations can be used to determine the distance 
modulus because either one of the PL relations would yield 
the same distance modulus after the appropriate metallic- 
ity correction. However, further study is needed to test 
these conclusions. 
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